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The past few years have seen a major advance in observational knowledge of high-energy (HE) pulsars. The Fermi Large 
Area Telescope (LAT) and AGILE have increased the number of known 7 -ray pulsars by an order of magnitude, its mem- 
bers being divided roughly equally among millisecond pulsars (MSPs), young radio-loud pulsars, and young radio-quiet 
pulsars. Many new and diverse emission characteristics are being measured, while radio and X-ray follow-up observations 
increase the pulsar detection rate and enrich our multiwavelength picture of these extreme sources. The wealth of new data 
has provided impetus for further development and improvement of existing theoretical pulsar models. Geometric light 
curve (LC) modelling has uncovered three broad classes into which HE pulsars fall: those where the radio profile leads, is 
aligned with, or lags the 7 -ray profile. For example, the original MSP and original black widow system are members of the 
second class, requiring co-located emission regions and thereby breaking with traditional notions of radio emission origin. 
These models imply narrow accelerator gaps in the outer magnetosphere, indicating copious pair production even in MSP 
magnetospheres that were previously thought to be pair-starved. The increased quality and variety of the LCs necessitate 
construction of ever more sophisticated models. We will review progress in global magnetosphere solutions which specify 
a finite conductivity on field lines above the stellar surface, filling the gap between the standard vacuum and force-free 
(FF; plasma-filled) models. The possibility of deriving phase-resolved spectra for the brightest pulsars, coupled with the 
fact that the HE pulsar population is sizable enough to allow sampling of various pulsar geometries, will enable much 
more stringent testing of future radiation models. Reproduction of the observed phase-resolved behavior of this disparate 
group will be one of the next frontiers in pulsar science, impacting on our understanding of particle acceleration, emis- 
sion, and magnetosphere geometry. One may now also study evolutionary trends of the measured or inferred quantities, 
and probe pulsar visibility and population properties such as radiation beam sizes of different pulsar classes, as well as 
the distribution of spin-down power, 7 -ray luminosity, conversion efficiency, spectral index, and cutoff energy across the 
population. Lastly, the recent detection of very-high-energy (VHE) pulsations from the Crab pulsar generated quite a few 
ideas to explain this emission, leading to an extension of standard models and possibly even a bridge between the physics 
of pulsars and pulsar wind nebulae (PWNe). 
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1 Introduction 

The incredible rate of 7 -ray pulsar detections by the Fenni 
Large Area Telescope (LAT) in the GeV band has brought 
tremendous progress to this field. The Second Pulsar Cat- 
alog (2PC; Abdo et al. 2013) details the properties of 117 
high-confidence detections above 100 MeV using 3 years 
of Fenni data, with 50 % of these pulsars having been dis- 
covered through 7 -ray pulsations or by radio searches of 
7 -ray sources. The population is distributed evenly among 
millisecond pulsars (MSPs), young radio-loud pulsars, and 
young radio-quiet pulsars. For a discussion of their multi- 
band profiles, off-peak emission properties, luminosities, 
spectra, and other population trends, see 2PC. Radio and X- 
ray follow-up observations have also played an important 
role in many of these detections, giving multiwavelength 
constraints on source energetics and viewing geometry. 

* Corresponding author: Christo. Venter@nwu.ac.za 


Many novel emission characteristics have come to light, 
stimulating a number of theoretical developments. We will 
summarize results from geometric light curve (LC) mod- 
elling (Sect. 2), focusing on both standard and new emission 
geometries and magnetospheric structures. Next, we will re- 
view results obtained from phase-resolved spectroscopy of 
the brightest 7 -ray pulsars (Sect. 3), and describe results of 
and prospects for TeV-domain pulsar science (Sect. 4). We 
lastly list some open questions in Sect. 5. 

2 Probing the magnetospheric and emi s sion 
geometry via light curve modelling 

The availability of high quality 7 -ray LCs has sparked a bar- 
rage of geometric LC modelling that attempts to inform us 
about the geometric properties of the magnetospheric emis- 
sion regions. This is a necessary first step that will guide 
the development of more physical radiation models. Below, 


Wiley Online Library 


© 2014 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim 


Astron. Nachr. /AN 335, No. 3 (2014) 


269 


we will first describe work using the standard models, for 
both young and old pulsars, before turning to newer emis- 
sion geometries and magnetosphere solutions. 

2.1 Standard emission geometries 

Some years after the solution by Deutsch (1955) of a mag- 
netic field rotating in vacuum (vacuum retarded dipole; 
VRD), Goldreich & Julian (1969) showed that the induced 
electric force vastly exceeds the gravitational pull on surface 
charges, so that the pulsar magnetosphere should be plasma- 
filled or “force-free” (FF; Spitkovsky 2006) rather than be- 
ing a vacuum. However, in a fully FF magnetosphere, the 
accelerating electric field (parallel to the magnetic field) 
is screened everywhere, and there can be no pulsed non- 
thermal emission. There should therefore be local devia- 
tions from this ideal assumption. This led to the develop- 
ment of a number of physical gap models. 

The polar cap (PC) model (Daugherty & Harding 1996) 
assumes charge injection from the PC, with radiation tak- 
ing place only at low altitudes (a few stellar radii) until 
secondary electron-positron pairs created from high-energy 
(HE) emission by relativistic primary particles screen the 
electric field at the so-called pair formation front. Taking 
the boundary condition of a zero electric field on the last 
open field line boundary properly into account leads to the 
pair formation front rising in altitude close to this boundary, 
forming a slot gap (SG) structure (Arons 1983; Muslimov & 
Harding 2003; Muslimov & Harding 2004a) in which parti- 
cles can be accelerated and emit radiation from the surface 
up to the light cylinder i?Lc (where the corotation speed 
equals the speed of light). The PC and SG models assume 
space-charge-limited flow, where charges are freely pulled 
from the stellar surface. The geometric two-pole caustic 
model (TPC; Dyks & Rudak 2003) resembles this SG struc- 
ture, with both magnetic poles being visible. The outer gap 
(OG) model (Cheng et al. 1986; Romani 1996) posits outer- 
magnetospheric gaps reaching from the null charge surface 
(NCS) where the Goldreich-Julian charge density is zero up 
to 7?lc- This represents a one-pole caustic model. ‘Caus- 
tic’ here refers to the bunching in phase of photons radiated 
on trailing magnetic field lines from a range of altitudes, 
due to the combination of time-of-flight delays, aberration, 
and magnetic field line sweepback (Morini 1983; Romani & 
Yadigaroglu 1995; Dyks & Rudak 2003). This is a feature of 
the standard outer-magnetospheric gap models, and is key to 
forming the 7-ray LC peaks. Additionally, when the voltage 
and magnetic field are too low to facilitate pair formation, 
particles may be accelerated in an unscreened open volume 
above the PC in the pair-starved polar cap (PSPC) model 
(Muslimov & Harding 2004b; Muslimov & Harding 2009). 

2.1.1 Younger pulsars 

Early on in the Fermi mission, Watters et al. (2009) pro- 
vided a practical visual summary of the beaming properties 


of 7-ray pulsars. Using PC, low-altitude TPC, and OG geo- 
metries, they provided maps indicating 7-ray peak multi- 
plicity and separation A as functions of the inclination and 
observer angles a and (j, and gap width w, and also separat- 
ing radio-loud and radio-quiet pulsars. Their “atlases” may 
be used in two ways: for a known pulsar geometry (a, C) 
from multiwavelength constraints, the 7-ray pulse proper- 
ties may be inferred for each of the different geometric mod- 
els. Conversely, given some observed HE pulse properties, 
the viewing geometry may be constrained, and one may 
sometimes discriminate between different geometric mod- 
els. Watters et al. (2009) noted that a plot of A vs. radio-to-7 
phase lag 5 (the “A-J relation”) may be a powerful discrim- 
inator between models, also constraining the radio emission 
properties. Lastly, they demonstrated the role of having ac- 
curate assessments of the correction factor /n, used when 
converting from the observed energy flux to intrinsic HE lu- 
minosity L^, as the 7-ray efficiency rj = L^/E (x fn (with 
E the spin-down power). 

Romani & Watters (2010) used Fermi LCs of 5 bright, 
younger pulsars, in addition to external multiwavelength 
constraints on the pulsar geometry (from X-ray pulsar wind 
nebula (PWN) torus fitting and radio polarization measure- 
ments), to constrain both the emission zone geometry and 
the magnetospheric structure. Using low-altitude TPC and 
OG geometries, as well as static, VRD, and pseudo-force- 
free (PFF) magnetic fields (the latter being constructed by 
adding a current-induced perturbation field to the VRD), 
they found that the OG model and PFF magnetosphere 
were statistically favored. Similarly, but following a differ- 
ent method, Li & Zhang (2011) used the VRD plus a per- 
turbed field (see Muslimov & Harding 2009) to study the 
effect of the latter on the LCs of the Vela pulsar, finding that 
a negative perturbation amplitude seems to be prefen'ed. 

Pierbattista et al. (2014) used PC, SG, and two OG geo- 
metries, in addition to conal and core radio models, in a 
VRD magnetosphere to model the LCs of 76 young and 
middle-aged LAT pulsars from 2PC. The bulk of the LCs 
favor outer-magnetospheric models (and high a and Q, al- 
though no single model best fits all LCs. PCs models favor 
low a and C^. Inclusion of radio LC fits as opposed to fitting 
7-ray profiles only leads to solutions closer to the diagonal 
on an a vs. plot, as required by radio visibility, so a and 
(j may otherwise be underestimated in the OG models. The 
values for /n are low for the PC model (given its narrow- 
beam geometry), /n ~ 1 for the SG, and /n ~ 0.5-0. 9 for 
the OG models (with radio-quiet pulsars generally requir- 
ing smaller values). Within errors, L.^ cx EE!'^ . There is less 
scatter for the modelled vs. E trend using the inferred 
/q than in 2PC (using /n = 1), so that this value may be 
an overestimation, especially for low-i5 objects. There is no 
evidence for evolution of a with age. Lastly, tentative in- 
verse trends for the 7-ray cutoff energy and spectral index 
vs. w were observed. 
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2.1.2 Millisecond pulsars 

Venter et al. (2009) geometrically modelled the LCs of 
the first 8 fermt-detected MSPs using standard 3D PSPC, 
TPC*, and OG models (and the radio cone model of Story 
et al. 2007). Since the bulk were best fit by TPC and OG 
models, copious pair production must be taking place to set 
up the narrow emission gaps, even in MSP magnetospheres 
that were previously thought to be pair-starved. Venter et al. 
(2009) also found /n ~ 1, similar to Watters et al. (2009). 
Venter et al. (2012) first noted that 7-ray pulsar LCs may 
be categorized into three broad classes (Class I to III): those 
where the radio leads the 7-ray profile (fit by standard OG 
/ TPC models), is phase-aligned with it, or lags the 7-ray 
peaks (exclusively fit by PSPC models). They developed 
two types of models for Class II pulsars, using co-located 
radio and 7-ray emission regions: altitude-limited TPC and 
OG (alTPC, alOG) models where the longitudinal gap ex- 
tent and position were varied given some restrictions, and 
low-altitude SG (laSG) models resembling hollow cones. 
LC fits favored the altitude-limited models, implying caus- 
tic origin of both the 7-ray and radio emission (Guillemot et 
al. 2012). The radio should be associated with depolariza- 
tion and rapid position angle swings (Dyks et al. 2004). 

Espinoza et al. (2013) introduced a different LC clas- 
sification scheme after studying the properties of 30 7-ray 
MSPs: (i) A-type MSPs exhibit phase-aligned radio and 7- 
ray peaks (equivalent to Class II MSPs; these also exhibit 
the largest magnetic fields at i?LC> steepest radio spectra, 
giant radio pulses in some cases, and low radio linear polar- 
ization which is possibly indicative of caustic radio emis- 
sion); (ii) N-type MSPs have non-aligned main radio and 
7-ray peaks, usually with single radio peaks; (iii) the main 
7-ray peak of W-type MSPs is out of phase with the main 
radio peak of a profile which consists of multiple peaks cov- 
ering a large phase range (“wide” peaks), possibly indicat- 
ing small a. The type N and W definitions intersect with the 
Class I and III designations discussed earlier. 

Johnson et al. (2014) systematically modelled the 
40 MSPs contained in 2PC using OG, TPC, laSG, and 
PSPC models, combined with a semi-empirical conal ra- 
dio model (Story et al. 2007), invoking core components 
when indicated by polarization measurements, and assum- 
ing a VRD field. They found that no single model best fits 
all LCs, although OG and TPC models provided best fits 
for the bulk of the MSP population (each best fits ~ 40 % 
of the MSPs). LCs with significant off-peak 7-ray emission 
are better fit by the TPC geometry, while LCs with little 
or no off-peak emission favour the OG geometry. Seven 
Class III MSPs are exclusively best fit by the PSPC model. 
This “mix” of results may point to some hybrid geometry 
that incorporates elements of the different standard models. 
The different LC classes pointed out by Venter et al. (2012) 
are not easily distinguished based on canonical pulsar vari- 

^ The “standard” or “conventional” TPC model assumes a maximum 
cylindrical radius of pcyl = 0.95/?lc for’ the emission region, while the 
“original” or “low-altitude” TPC assumes pcy\ = 0.75i?LC- 


ables (such as spin period, magnetic field, or E) alone, but 
may rather reflect more complex electrodynamical environ- 
ments for the different pulsars. The best-fit a occurred over 
a broad range of values, which may be a result of their spin- 
up evolution and may also reflect the relatively wider ra- 
dio beams and smaller magnetospheres of MSPs, making 
them visible over a larger region of phase space than their 
younger counterparts. This implies that there are probably 
not many radio-quiet MSPs. A clustering at large ( is ob- 
served, ascribed to the bright caustic emission lying near 
the spin equator which is usually sampled to form the HE 
peaks. A preponderance of ~ 90° is expected if pul- 
sar spin axes are distributed randomly with respect to the 
Earth line of sight. OG models prefer larger a, since these 
gaps may be invisible for small a (depending on w). Eor 
MSPs, Lj oc E, confirming the expectation that older pul- 
sars’ magnetospheres are more pair-starved, enabling them 
to reach relatively larger 7 ~ 10 %. 

2.2 Other magnetospheres or emission geometries 
2.2.1 TPC / OG models in an FF magnetosphere 

While numerous geometric LC studies have been done 
using the VRD as the basic magnetospheric structure 
(Sect. 2.1), the FF magnetic structure should be more re- 
alistic in that it takes into account the effect of charges and 
currents. The FF magnetosphere satisfies pE+j xB jc = 0, 
with p and j the charge and current densities, such that 
E ■ B = 0, implying that particle acceleration or emission 
is precluded (see, e.g., Contopoulos et al. 1999; Gruzi- 
nov 2005; Komissarov 2006; McKinney 2006; Spitkovsky 
2006; Timokhin 2006; Kalapotharakos et al. 2012a). 

Bai & Spitkovsky (2010) presented model LCs using 
the conventional TPC and OG models, but in an FF mag- 
netosphere. They assumed constant emissivity along parti- 
cle trajectories in the lab frame and found that the FF TPC 
model cannot produce sharp double-peaked LCs due to the 
increased PC size compared to that of the VRD. This sup- 
presses caustic formation close to the stellar surface, so that 
all caustics are produced close to i?LC- For the OG case, 
emission is gathered only from open field lines that cross 
the NCS, and this generally eliminates one peak from the 
resulting LC. Next, Bai & Spitkovsky (2010) investigated 
a “separatrix layer model” in the FF magnetosphere (very 
similar to the FF TPC, but slightly interior to it) where 7- 
ray emission originates in a thin layer in the open-field-line 
region just inside the separatrix (current sheet) which sepa- 
rates open and closed field lines. Field lines with footpoints 
close but interior to the PC rim are involved, and emis- 
sion extends from the stellar surface to beyond i?LC ■ A new 
caustic (one from each pole) appears which is due to “sky 
map stagnation” (SMS), where emission originating close 
to 7 ?lc arrives at the same observer phase. Double-peaked 
LCs (mostly at large a and Q are a generic feature due 
to caustics from both poles that form in the outer magne- 
tosphere, where the field lines approach a split-monopolar 
geometry. The SMS effect becomes more pronounced as 
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one chooses radiating field lines closer to the magnetic 
axis.^ 

Harding et al. (2011) modelled TPC LCs in VRD and 
FF magnetospheres, investigating the effect of an azimuthal 
modulation in emissivity, since the PCs in these magneto- 
spheres are offset from the magnetic axis, inducing an azi- 
muthally asymmetric electric field (Harding & Muslimov 
2011). LCs incorporating this modulation exhibit much 
lower off-peak emission. The FF LC peaks are shifted later 
in phase (vs. VRD LCs), since the FF PCs are larger and the 
field lines are more swept back. They found that the Vela 
LCs favor the asymmetric VRD LCs, possibly indicating 
that the real pulsar magnetosphere solution is closer to the 
VRD than the FF geometry. DeCesar et al. (2014) obtained 
best-fit TPC and OG LCs for four bright pulsars (Vela, Crab, 
Geminga, and PSR J0007-(-7303), using VDR and FF fields. 
No single emission model or magnetic geometry was pre- 
ferred, pointing to yet more complex emission geometries. 

2.2.2 Annular gap model 

Qiao et al. (2004) proposed an annular gap model to explain 
pulsar emission. On the stellar surface, an annulus is formed 
by the footpoints of the last open field lines (PC rim) and 
those of the critical field lines which intersect with the NCS 
at 7?lc- This bounded region is referred to as the annular 
gap (AG), while the region interior to the AG is called the 
core gap (CG). The annulus shape is a strong function of a 
(assuming a dipolar magnetic field), and may be a source of 
relativistic particles that radiate between the AG and NCS. 

Du et al. (2010) presented geometric HE LC modelling 
(no radio modelling) using a single-pole, intermediate- 
height 3D AG (favoring short period-pulsars). They as- 
sumed two Gaussian distributions for the emissivity: one 
along the field lines (peaking near the NCS), and one across 
the gap in the latitudinal direction (peaking near the AG 
center). By including aberration and retardation and also 
taking (, constraints from PWN torus fitting (Ng & Romani 
2008) into account, they obtained reasonable fits for the LCs 
of 6 Fermi pulsars, noting that their model emission from 
younger pulsars peaked near the NCS, while it originated 
from much lower altitudes for MSPs (implying smaller 
beam solid angles AH). For Vela, Du et al. (2011) invoked 
HE emission from the AG to fit the two main peaks (PI 
and P2), while the third “bump” (P3) was fit using the CG. 
They fit multiband LCs by adjusting the dominant emission 
height for each energy band. Radio emission was assumed 
to come from a narrow region near i?LC- Using a 2D global 
electric field, they calculated phase-resolved spectra assum- 
ing synchro-curvature radiation (synchro-CR) from primary 
particles and synchrotron radiation (SR) from secondaries, 
assuming power-law spectra and adjusting AH. There are 
many free parameters in this approach. 

^ See also Contopoulos & Kalapotharakos (2010) who studied signa- 
tures in LCs and polarization profiles of CR originating in the strong cur- 
rent sheet that develops at the tip of the pulsar closed line region using 
time-dependent 3D FF magnetosphere. 


Du et al. (2012) next modelled the radio-to-TeV emis- 
sion from the Crab pulsar in seven phase bands (plus phase- 
averaged) using the single-pole AG model. They used a ID 
accelerating potential along the field lines, as well as power- 
law distributed primaries and two populations of pairs, in- 
duced by CR and inverse Compton scattering (ICS) respec- 
tively. SR from pairs dominated the X-ray to soft 7-ray 
band; CR and SR from primaries the 7-ray band; and ICS 
from pairs (on thermal surface and SR photons) the TeV 
band. In order to constrain the emission geometry and loca- 
tions, they used prior geometric LC fitting invoking the AG 
to model the peaks, and the CG to model the bridge emis- 
sion. Du et al. (2013) modelled the radio and 7-ray LCs 
of three MSPs (PSR 10034-0534, PSR 10101-6422, and 
PSR 10437-4715) using a single-pole AG plus CG model. 
The challenging three-peaked 7-ray LC of PSR 10101-6422 
could be fit invoking emission from both the AG and two 
zones in the CG, coming from different emission altitudes. 

2.2.3 Striped wind model 

Building on previous work (Coroniti 1990; Michel 1994), 
Petri (2009) considered the striped wind as an alternative 
emission region outside i?LC to generate HE emission. He 
used an asymptotic oblique split monopole solution (Bo- 
govalov 1999), noting that relativistic beaming effects give 
rise to pulsed emission from this region. Calculating ICS 
from an electron-positron pair wind on cosmic microwave 
background photons, Petri could reproduce the LCs and 
phase-resolved spectra of Geminga using EGRET data. 
Petri (201 1) presented calculations of LCs for Fermi pulsars 
(neglecting energy dependence), also including radio emis- 
sion from the PCs, modeled as a Gaussian emissivity profile 
peaking on the magnetic axis. Petri noted the analytic rela- 
tion cos(7T A) = I cot C cot a\, first found by Kirk (2005), as 
well as an expression for the A-J relation (5 rj 0.5 — 0.5A), 
interpreting the phase lag 5 as a time-of-fiight effect. Two 
regions of roughly constant /o(«, C) ~ {0-4, 1.9} were ob- 
tained. The S in this model are generally larger than seen 
in the data. Another application considers 7-ray SR from a 
relativistically hot current sheet in the striped wind, includ- 
ing some magnetic reconnection (Petri 2012). 

2.2.4 Dissipative magnetospheres 

Newer global magnetosphere solutions specify a (constant) 
finite conductivity a on field lines above the stellar sur- 
face {E ■ B ^ 0 so that there may be particle acceleration), 
thereby filling the gap between the vacuum and FF mod- 
els; see Kalapotharakos et al. (2012b) and Li et al. (2012). 
Kalapotharakos et al. (2012b) considered an orthogonal ro- 
tator, and followed two approaches to calculate HE LCs: the 
usual geometric approach (Sect. 2.1), and a particle trajec- 
tory approach (since in these magnetospheres, particles do 
not follow the magnetic field lines in the corotating frame) 
similar to that used in the separatrix model (Sect. 2.2.1). 
For larger a (larger field sweepback), there is a progressive 
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LC shift to larger phases and peak broadening. For the sec- 
ond approach, narrowest pulses are found for high a, and 5 
generally decreases with for lower a values. When exclu- 
sively considering emission close to f?LC in these magne- 
tospheres, improved predictions for the A-5 trend may be 
obtained (Kalapotharakos et al. 2014). 

3 Phase-resolved spectroscopy 

The vast number of accumulated statistics from the bright- 
est 7-ray pulsars is now allowing us to probe the emission 
environment and processes on a much deeper level than 
just using LC modelling. For example, Abdo et al. (2010) 
presented phase-resolved spectra of the Vela pulsar, show- 
ing the variation of photon index and cutoff energy Ec as 
a function of phase (j). Although this has been done previ- 
ously (e.g., Kanbach et al. 1994; Fierro et al. 1998), much 
finer phase bins may now be attained and more complicated 
spectral shapes may be fit to data in a particular phase bin. 

DeCesar et al. (2014) studied four bright Fermi pulsars 
using 30 months of data. The LC peaks are thought to be 
due to caustic emission (Sect. 2.1), with photons originat- 
ing from very different emission radii being contained in 
any particular phase bin. Therefore, the inferred reflects 
the combination of emitted spectra from various positions 
in the magnetosphere (sampling different local conditions, 
e.g., magnetic field and curvature radius pc). This compli- 
cates the disentanglement of altitudinal behavior of the ac- 
celerating electric field fJ| | (r) when one only finds an indi- 
cation of the combined effect of many different local i?| | and 
Pc contributing to the spectral shape at a particular phase (j) 
through the measured Ec{4>)- Nevertheless, when assuming 
that CR is the dominant emission process in the GeV band 
and that the CR reaction limit is valid, one may constrain 
the implied i5|| (0) using the measured Ec{4>) and modelled 
^max(^), ^ particular inferred pulsar geo- 

metry (from the LC modelling) and magnetosphere, and fi- 
nally infer E\ | (rniin) using the modelled minimum emission 
radius r-min (<)>)■ This in turn constrains the magnetospheric 
structure, implying that high-altitude field lines should be 
straighter (higher pd compared to the vacuum case if E\ \ is 
not to exceed the local magnetic field (as it should be, since 
it is induced by this magnetic field in the first place). The 
grand challenge for future physical emission models will be 
to reproduce both the energy-dependent LCs as well as the 
phase-dependent spectral behavior seen in a population of 
bright 7-ray pulsars, each having unique geometries, as this 
will yield the most stringent constraints to date on particle 
acceleration, emission, and magnetic field geometry in the 
pulsar magnetosphere. 

4 Very-high-energy pulsar science 

Many of the Fermi-detected pulsars exhibit exponentially 
cutoff spectra (2PC). This has been interpreted as CR by 


relativistic particles moving along curved field lines. Typ- 
ical pulsar parameters result in cutoffs around a few GeV, 
exactly what Fenni has seen (e.g.. Venter & de Jager 2010). 
While earlier models predicted that these cutoffs may reach 
100 GeV (Bulik et al. 2000; Harding et al. 2002) in MSPs, 
later PSPC MSP models predicted lower cutoffs ~ 10- 
50 GeV (Harding et al. 2005; Venter & de Jager 2005; 
Fr^ckowiak & Rudak 2005). The SG model applied to the 
Crab pulsar predicted cutoffs up to a few GeV (Harding et 
al. 2008), while the OG model could explain 7-rays of up to 
25 GeV (e.g., Hirotani 2008; Tang et al. 2008). Given this 
theoretical context, it was very surprising to detect pulsed 
emission at tens to hundreds of GeV from the Crab pulsar 
(Aliu et al. 2008,Saito 2010, Aleksif et al. 201 1, Aliu et al. 
201 1, Aleksi§ et al. 2012), indicating a necessary extension 
of the standard pulsar models. 

Spectral cutoffs are determined by the maximal particle 
acceleration, as well as (one- and two-photon) absorption 
of energetic 7-rays (e.g., Aliu et al. 2008). One needs an 
unrealistically high value for the accelerating electric field 
(Fj I ^ B) and/or radius of curvature (pc ~ ^?lc) to explain 
the observed VHE emission when invoking CR with a spec- 
tral cutoff or break around 100 GeV (Aliu et al. 2011). 

A revised OG model by Hirotani may produce inverse 
Compton (IC) radiation of up to 400 GeV by secondary 
and tertiary pairs on infrared to ultraviolet (UV) photons to 
explain the observed Crab spectrum (Aleksij et al. 2012). 
Extending the calculations within an SG model, assuming 
non-isotropic IC of SR photons by the pair plasma, may 
prove useful. This approach is however very dependent on 
the assumed electrodynamics and magnetic field geometry, 
making the spectral prediction somewhat uncertain. Lyu- 
tikov et al. (2012) proposed a synchrotron self-Compton 
model to explain the VHE emission. They assume that the 
electron-positron plasma upscatter UV SR photons in the 
Klein-Nishina regime and obtain a reasonable reproduction 
of the pulsed VHE spectrum. They also argue for a sta- 
tistical correlation between X-ray and 7-ray photons on a 
short timescale. The difference in HE and VHE LCs, and 
the flattening of the spectrum with energy motivated Aharo- 
nian et al. (2012) to invoke a low-energy, mono-energetic, 
ultrarelativistic wind dominated by kinetic particle energy 
that is very rapidly accelerated close to i?LC to upscatter 
(mainly X-ray) photons from the pulsar surface and/or mag- 
netosphere via anisotropic IC scattering. Other possibilities 
beyond extensions to the standard scenarios include the AG 
model by Du et al. (2012), or SR by primary electrons with 
very high Lorentz factors (7 ~ 10®; Chkheidze et al. 201 1). 

5 Conclusion 

With the availability of high-quality 7-ray data the empha- 
sis has now shifted to model refinement, entailing a move 
beyond the standard vacuum-based (or even EF) magne- 
tosphere models toward self-consistent, dissipative magne- 
tospheres that determine both the geometry and position- 
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dependent emissivity of the emitting region from first prin- 
ciples. The question of choosing a best model should there- 
fore change from a “TPC vs. OG argument” to constraints 
on the plasma properties (e.g., density and conductivity), 
also taking into account the effect of pulsar evolution. We 
are thus moving toward constraining the actual pulsar elec- 
trodynamics, and not just geometry. However, integrating a 
physics-motivated prescription for the conductivity (the mi- 
crophysics of pair cascades) into the dissipative magneto- 
sphere models is a necessary next step. Lastly, a population 
synthesis approach (e.g., Watters & Romani 2011; Takata et 
al. 2011; Pierbattista et al. 2012) which includes more real- 
istic magnetospheres and emission geometries should be a 
key avenue for future progress. 

We conclude this review by listing some open questions: 

1. There is a rich variety of LCs in the newly-detected 
7 -ray pulsar population; is this solely due to geometry, 
or may distinct magnetospheric conditions and evolu- 
tion also play a role? 

2. There is no overall best-fit geometrical LC model; this 
points to some hybrid or new emission geometry. What 
would such a geometry be? 

3 . Preliminary trends are seen between geometric / spectral 
vs. pulsar parameters. Will these trends strengthen with 
more pulsar discoveries, or do we need to study a larger 
phase space, including noncanonical pulsar parameters? 

4. Which is the preferred magnetospheric structure? What 
is the effect of the current sheet? Does magnetic recon- 
nection play a role to dissipate particle energy and create 
emission? 

5. The behavior of spectral parameters with phase is highly 
complex; how is this explained? Will similar behavior 
be seen when phase-resolved spectroscopy of dimmer 
pulsars becomes viable? 

6. What is the origin of pulsed TeV photons from pulsars? 

Acknowledgements. CV is supported by the South African Na- 
tional Research Foundation. AKH acknowledges support from the 
NASA Astrophysics Theory Program. CV and AKH acknowledge 
support from the Fermi Guest Investigator Program. 

References 

Abdo, A.A., et al. 2010, ApJ, 713, 154 
Abdo, A.A., et al. 2013, ApJS, 208, 17 (2PC) 

Aharonian, F.A., et al. 2012, Nature, 482, 507 

Aleksi 9 , J., et al. 2011, ApJ, 742, 43 

Aleksif, J., et al. 2012, A&A, 540, A69 

Aliu, E., et al. 2008, Science, 322, 1221 

Aliu, E., et al. 2011, Science, 334, 69 

Arons, J. 1983, ApJ, 266, 215 

Bai, X.-N., & Spitkovsky, A. 2010, ApJ, 715, 1282 

Bogovalov, S.V. 1999, A&A, 349, 1017 

Bulik, T, Rudak, B., & Dyks, J. 2000, MNRAS, 317, 97 

Cheng, K.S., Ho, C., & Ruderman, M. 1986, ApJ, 300, 500 

Chkheidze, N., et al. 2011, ApJ, 730, 62 

Contopoulos, I., Kazanas, D., & Eendt, C. 1999, ApJ, 511, 351 

Contopoulos, I., & Kalapotharakos, C. 2010, MNRAS, 404, 767 

Coroniti, E.V. 1990, ApJ, 349, 538 


Daugherty, J.K., & Harding, A.K. 1996, ApJ, 458, 278 

DeCesar, M.E., et al. 2014, in preparation 

Deutsch, A.J. 1955, Ann. d’Astrophys., 18, 1 

Du, Y.J., et al. 2010, MNRAS, 406, 2671 

Du, Y.J., Han, J.L., Qiao, G.J., & Chou, C.K. 2011, ApJ, 731, 2 

Du, Y.J., Qiao, G.J., & Wang, W. 2012, ApJ, 748, 84 

Du, Y.J., Qiao, G.J., & Chen, D. 2013, ApJ, 763, 29 

Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 

Dyks J., Harding, A.K., & Rudak, B. 2004, ApJ, 606, 1 125 

Espinoza, C.M., et al. 2013, MNRAS, 430, 571 

Eierro, J.M., et al. 1998, ApJ, 494, 734 

Erjckowiak, M., & Rudak, B. 2005, Adv. Space Res., 35, 1152 

Goldreich, P, & Julian, W.H. 1969, ApJ, 157, 869 

Gruzinov, A. 2005, Phys. Rev. Lett., 94, 021 101 

Guillemot, L., et al. 2012, ApJ, 744, 33 

Harding, A.K., Muslimov, A.G., & Zhang, B. 2002, ApJ, 576, 366 
Harding, A.K., Usov, V.V., & Muslimov, A.G. 2005, ApJ, 622, 531 
Harding, A.K., Stern, J.V., Dyks, J., & Erjckowiak, M. 2008, ApJ, 
680, 1378 

Harding, A.K., DeCesar, M.E., Miller, M.C., Kalapotharakos, C., 
& Contopoulos, I. 201 1 Fermi Symp. Proc., eConf Cl 10509 
Harding, A.K., & Muslimov, A.G. 2011, ApJ, 743, 181 
Hirotani, K. 2008, ApJ, 688, L25 
Johnson, T.J., et al. 2014, in preparation 

Kalapotharakos, C., Contopoulos, J., & Kazanas, D. 2012a, MN- 
RAS, 420, 2793 

Kalapotharakos, C., et al. 2012b, ApJ, 754, LI 
Kalapotharakos, C., et al. 2014, in preparation 
Kanhach, G., et al. 1994, A&A, 289, 855 
Kirk, J.G. 2005, Mem. Soc. Astron. Ital., 76, 494 
Komissarov, S.S. 2006, MNRAS, 367, 19 
Li, X., & Zhang, L. 2011, ApJ, 743, 113 
Li, J., Spitkovsky, A., & Tchekhovskoy, A. 2012, ApJ, 746, 60 
Lyutikov, M., Otte, N., & McCann, A. 2012, ApJ, 754, 33 
McKinney, J.C. 2006, MNRAS, 368, L30 
Michel, F.C. 1994, ApJ, 431, 397 
Morini, M. 1983, MNRAS, 202, 495 
Muslimov, A.G., & Harding, A.K. 2003, ApJ, 588, 430 
Muslimov, A.G., & Harding, A.K. 2004a, ApJ, 606, 1 143 
Muslimov, A.G., & Harding, A.K. 2004b, ApJ, 617, 471 
Muslimov, A.G., & Harding, A.K. 2009, ApJ, 692, 140 
Ng, C.-Y, & Romani, R.W. 2008, ApJ, 673, 411 
Petri, J. 2009, A&A, 503, 13 
Petri, J. 2011, MNRAS, 412, 1870 
Petri, J. 2012, MNRAS, 424, 2023 
Pierbattista, M., et al. 2012, A&A, 545, A42 
Pierbattista, M., et al. 2014, ApJ, in preparation 
Qiao, G.J., et al. 2004, ApJ, 606, L49 
Romani, R.W. 1996, ApJ, 470, 469 
Romani, R.W., & Watters, K.P 2010, ApJ, 714, 810 
Romani, R.W., & Yadigaroglu, I.-A. 1995, ApJ, 438, 314 
Saito, T. 2010, PhD thesis, Ludwig-Maximilians-Universitat 
Miinchen, Germany 
Spitkovsky, A. 2006, ApJ, 648, L5 1 

Story, S.A., Gonthier, P.L., & Harding, A.K. 2007, ApJ, 671, 713 

Takata, J., Wang, Y, & Cheng, K.S. 2011, ApJ, 726, 44 

Tang, A.P.S., et al. 2008, ApJ, 676, 562 

Timokhin, A.N. 2006, MNRAS, 368, 1055 

Venter, C., & de Jager, O.C. 2005, ApJ, 619, L167 

Venter, C., Harding, A.K., & Guillemot, L. 2009, ApJ, 707, 800 

Venter, C., & de Jager, O.C. 2010, ApJ, 725, 1903 

Venter, C., Johnson, T.J., & Harding, A.K. 2012, ApJ, 744, 34 

Watters, K.P., et al. 2009, ApJ, 695, 1289 

Watters, K.P., & Romani, R.W. 2011, ApJ, 727, 123 


www.an-joumal.org 


© 2014 WILEY- VCH Verlag GmbH & Co. KGaA, Weinheim 


